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ABSTRACT 



Motivated by Papadakis (2005a, 2005b), we study a Chermnykh-like problem, 
in which an additional gravitational potential from the belt is included. In 
addition to the usual five equilibrium points (three collinear and two triangular 
points), there are some new equilibrium points for this system. We studied the 
conditions for the existence of these new equilibrium points both analytically 
and numerically. 
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1. Introduction 

The interest in the Chermnykh's problem (Chermnykh 1987) has been renewed by 
K.E. Papadakis recently. For example, Papadakis (2004) studied the symmetric motions 
near the three collinear equilibrium points in three dimensional space. Papadakis (2005a) 
numerically investigated the equilibrium points and the zero-velocity curves under the 
assumption that the mass parameter is fixed but the angular velocity parameter is allowed 
to vary continuously. The stability of the periodic orbits are determined and the network of 
the orbital families is explored. 

Moreover, Papadakis (2005b) numerically studied the asymmetric periodic orbits near 
the triangular equilibrium points under the assumption that the angular velocity varies 
and for the Sun- Jupiter mass distribution. Remarkably, that paper provides the analytic 
determination of the initial conditions of the long- and short-period Trojan families around 
the equilibrium points. 

On the other hand, astronomers claimed that it is common to have circumbinary discs 
for binary systems and asteroid belts for planetary systems. Some studies showed that these 
belt-like structures shall influence the mathematical properties of these dynamical systems. 
The number and locations of equilibrium points and also topology of solution curves might 
become very different when the influence from the belt is considered. For example, Jiang 
& Yeh (2003, 2004) considered the influence from the belt for planetary systems and found 
that the probability to have equilibrium points around the inner part of the belt is larger 
than the one near the outer part. Their results can be used to explain the observational 
configuration of Kuiper Belt Objects of the outer Solar System. 

Jiang & Yeh (2006) have studied a Chermnykh-like problem in which the mass 
parameter fi is set to be 0.5. We now study the systems with any possible mass-ratio of 
central binaries both analytically and numerically. We will focus on conditions with which 
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the new equilibrium points can exist for two different density models of the belts. We can 
then understand the locations of these new equilibrium points for given parameters. We 
present our models in Section 2, the results for the existence of new equilibrium points are 
in Section 3. Section 4 is about the stability analysis and Section 5 concludes the paper. 



The Models 



We consider the motion of a test particle influenced by the gravitational force from the 
central binary and the circumbinary belt. We assume that two masses of the central binary 
are mi and m 2 and choose the unit of mass to make G{m\ + m 2 ) = 1. If we define that 



m 2 



H = 

m 1 + m 2 

then the two masses are /ii = Gm\ — 1 — fj, and /i 2 = Gm 2 = /i. The separation of central 
binary objects is set to be unity. In this paper, we assume that /i < 0.5. The location of 
two masses are always at (1 — /i, 0, 0) and (— /i, 0, 0). 

The equation of motion is (Jiang & Yeh 2006) 



du 
dt 

dv 
dt 



ox ox 



(1) 



-2nu - 

oy 



dV 
dy ' 



where n is the central binary's angular velocity (Please note that we keep n as a parameter 
from here until Property 3.2 but set it to be 1 as a simplification starting from Eq.(20), i.e, 
Section 3.1 Model A.), 

n 2 _ . 1 - ,/ 

(2) 



£/* = -^ 2 + y 2 )-^ 

Z T\ 



r 2 



(3) 
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(4) 



and V is the potential from the belt. We assume that V is radial symmetric, so V depends 
on the radial distance r, where r = \Jx 2 + y 2 . Hence, 



av 

dx 



■ dV 
dr 



dV _ ydV 
dy r dr ' 



(5) 



We substitute Eq. (2) and Eq. (5) into Eq.(l) and have the following system: 



dx 
dt 



U 



dt u 

4j± = 2nv + n 2 x 



(1— jl)(x+jf) fl(x + fl—l) 



x dV 
r dr 



(6) 



f = -2nu + n 2 y - 



y(i-M) 



r 2 



ydV 
r dr ' 



where r\ and r 2 are defined in Eqs.(3)-(4). We can consider any density profile of the belt 
by giving a formula of V in Eq. (6) as long as V is radial symmetric. In this paper, we 
study two different models of the belt. 



2.1. Model A: Miyamoto-Nagai Profile 

In Model A, we introduce the belt potential, from Miyamoto & Nagai (1975), which is 

V(r, z) = - ^^^ => (7) 

\jr 2 + (a + Vz 2 + b 2 ) 2 

where M b is the total mass of the belt and r 2 = x 2 + y 2 . In this formula, a and b are 
parameters which determine the density profile of the belt. The parameter a controls the 
flatness of the profile and can be called "flatness parameter" . The parameter b controls the 
size of the core of the density profile and can be called "core parameter" . When a = b = 0, 
the potential equals to the one by a point mass. In general, the density distribution 
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corresponding to the above V(r,z) in Eq.(7) is, as in Miyamoto & Nagai (1975), 

- ar ' 2 + ( a + 3 ^ 2 + &2 )( a + y^ 2 + ft 2 ) 2 

If we define T = a + 6, from Eq. (7) we have 



(8) 



dV 



r,0) 



M b r 



M b r 



(9) 



dr v '-' (r 2 + (a + 6) 2 ) 3 / 2 ( r 2 + T 2 ) 3 / 2 ' 
where we set z = since we only consider the orbits on the x — y plane. We substitute Eq. 
(9) into Eq. (6) and have the following system: 



u 



dx 
dt 

^ = 2nv + n 2 x - 
§ = -2nu + n 2 y 



(l-fj,)(x+fj,) n(x+n—l) 



M b x 

( r 2 +T 2)3/2 



(10) 



M b y 

( r 2 +T 2)3/2 • 



We can see that neither the core parameter or the flatness parameter appear in the 
equations of motion. The dynamics of the system only depends on the summation of a and 
b, i.e. T. 



2.2. Model B: Power-Law Profile 



The belt is a annulus with inner radius r\ and outer radius r a , where and r Q are 
assumed to be constants. The density profile of the belt is 

when r < r, or r > r a , 



p{r) 



cos 



c 



f (ri-ri) } } when Ti<r <r u 
when ri < r < r 2 , 



(11) 



k ^{cos 



tt (r-r 2 ) 
2 (r D -r 2 ) 



1 2 

J- when 



r 2 < r < r oj 



where r = y/x 2 + y 2 , c is a constant completely determined by the total mass of the belt 
and we set p = 2 (Lizano & Shu 1989), r\ = r\ + 0.1, r 2 = n + 0.8, and r Q = r 2 + 0.1 for all 
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numerical results. Hence, for p — 2, the total mass of the belt is 



i-2-k pro i rri 1 

M b = p(r')r'dr'd(f) = 2nc { - 

JO JTi JTi T 



COS 



7r (r — r\) 
2 (n -ri) 



dr' + ln(r 2 /ri) 



r Q 1 



7r (r - r 2 ) 



dr' 



2 (r D - r 2 ) 

The gravitational force f b from the belt is (Jiang & Yeh 2004) 



(12) 



E(0 + ^(0 



dr', 



where £ = 2y/rr' '/(r + r'), F and E are elliptic integrals of the first kind and the second 
kind. We substitute Eq. (13) into Eq.(6) and have the following system: 



(13) 



dx 
dt 

dy 
~dt 



U 



^-2nv + n 2 x - - 



(if 
cii.' 



2x 
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f = -2rax + n 2 y 



%Sr r °p(r'y 



E _|_ F 



dr' 



E _|_ F 



dr'. 



(14) 



3. The New Equilibrium Points 



It is well-known that there are five equilibrium points, i.e. three collinear and two 
triangular points, for the classical restricted three-body problem. When there are more than 
five equilibrium points for the system we study here, we claim that the new equilibrium 
points exist. 



In general, for System (6), equilibrium points (x e ,y e ) satisfies f(x e ,y e 
g(x e ,y e ) = 0, where 

(l — fi)(x + n) /jl(x + ijl — 1) xdV 



and 



f(x,y) = n 2 x - 
g(x,y) = n 2 y - 



rf 

2 y(i-^) w v dv 



rf 



r dr ' 



rf 



(15) 
(16) 
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For convenience, we define 

h(y) = n 2 



+ y* 



-3/2 ^ l QV 

r dr 



and 



k(x) = n 2 x — 



(l-(i)(x + n) fi{x + fi-l) xdV 



\x + n\ 



\x + fi — 1| 3 r dr 



(x,0) 



(17) 



(18) 



We have the following properties: 



Property 3.1 

If (x e ,y e ) is an equilibrium point of System (6), then we have 
either (1) x e = 1/2 — /i and y e satisfies h(y) = 
or (2) x e satisfies k(x) =0 and y e = 0. 

Proof: Suppose that (x e ,y e ) is an equilibrium point, thus it satisfies f(x,y) = and 
g(x,y) = 0. From Eq.(16), we have 

I — fj> n 1 dV 



n 



rf r| v dr 



= 0. 



- Oe,2/e) 



Hence, y e = or n 



2 _ i-/x _ u _ i ay 



0. We discuss these two cases separately: 



(I)n 



2 _ i-M _ m _ i av 

r? r? r 9r 



0: 



Since (x e ,y e ) is an equilibrium point, that is f(x e ,y e ) = 0, we have 

1 — ji fx 1 dV 



= f(x e ,y e ) = x e 



n 



rf r dr 



rf 



rl 



1 1 



(19) 



Thus, ri = r 2 , i.e. (x e + /i) 2 + y 2 = (x e + fi — l) 2 + y 2 . Hence x e = 1/2-/1. We have 

-3/2 

( = n'- 7 + 3T 

' 1 '2 



ri = r 2 = Jl/A + y 2 e and thus, n 2 - ^ - & = n 2 - 



Therefore, x e = 1/2 — fi 



and y e satisfies fo(y) = for the equilibrium point (x e ,y e ). 



(II) y e = 0: 

f(x e ,y e ) = f(x e ,0) = k(x e ) = for the equilibrium point (x e ,y e ). Thus, x e satisfies 
k(x) = and y e = 0. □ 
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Property 3.2 

(A) If y e satisfies h(y) =0, then (1/2 — /i, y e ) is the equilibrium point of System (6). 

(B) If x e satisfies k(x) = 0, then (x e , 0) is the equilibrium point of System (6). 

Proof of (A): Suppose that y e is one of the roots of h(y), i.e. h(y e ) = and we set 
x e = l/2-fi. Because #(1/2 - fi, y e ) = y e h(y e ) = and f(l/2 - fi, y e ) = (| - fij h(y e ) = 0, 
(x e ,y e ) is the equilibrium point of System (6). □ 

Proof of (B): Suppose that x e is one of the roots of k(x), i.e. k(x e ) = and we set 
y e = 0. Since y e = 0, it is trivial that g(x e , y e ) = 0. Because k(x e ) = 0, f(x e , 0) = k(x e ) = 0. 
Thus, (x e ,y e ) is the equilibrium point of System (6). □ 

Because of the above two properties, in stead of searching the roots for two variable 
functions f(x, y) and g(x, y) to determine all equilibrium points on the x — y plane, we only 
need to find the roots of one variable functions h(y) and k(x) to get all the equilibrium 
points of System (6). 



3.1. Model A 



We now study the equilibrium points when V is Miyamoto-Nagai profile. From 
Property 3.1, equilibrium points should be either on the a>axis or on the line: x = 1/2 — /i 
of x — y plane. As a simplification, please note that we set n — 1 hereafter. For equilibrium 
points on the x-axis, i.e. (x e , 0), we have 

(1 - fj)(x e + jj) fx(x e + fx - i) x e M b 



k^Xgj x e 



0. 



\x e + /x| 3 \x e + fi — 1| 3 {x\ + T 2 ) 3 / 2 

On the other hand, for those equilibrium points on the line: x = 1/2 — /i (i.e. 

[1/2 — fi, y e )), and thus r = Ml/2 — fi) 2 + y 2 , we have 



(20) 



%e) = 1 



Ve 



-3/2 



1 



fi 



y 2 e + T 2 



-3/2 



= 0. 



(21) 
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Property 3.3 

There is one and only one yi > such that h(y~i) = 0, and only one y 2 < such that 
h{y~2) = 0. That is, there are two equilibrium points on the line: x — 1/2 — /i for System 
(10). 



Proof: We define 



PM = 



- 4 +y 



-3/2 



and Q 1 (y) = -M b 
so from Eq.(21), we have 



-3/2 





r 1 2 


-3/2 




-3/2 


h(y) = 1 - 


A +y _ 


-M b 




= -P 1 (y) + Q 1 (y)=0 



(22) 
(23) 

(24) 



From Eqs. (22)-(23), we have 
Q' 1 (y)=3yM b ' ' 



+y z + T' 



-5/2 



and P[{y) = -3y 



-5/2 



Since P[(y) < for y > 0, -Pi(y) is a monotonically decreasing and continuous function for 
any y G (0, oo). Q[(y) > for y > 0, thus Qi(y) is a monotonically increasing and continuous 
function for any y G (0, oo). Because -Pi(O) > > <5i(0) and lim^oo P\{y) < lim^oo Qi(y), 
there exists a unique point yi > such that Pi(y~i) = Qi(y~i), that is h(y~i) = 0. 

By the similar method, we find that there is a unique point y~2 < such that ^(2/2) = 0. 
Hence, there are two equilibrium points on the line: x — 1/2 — /i. □ 

In the classical restricted three-body problem, there are three collinear points and two 
triangular points, i.e. the points on the line: x — 1/2 — /i. Therefore, Property 3.3 shows 
that there is no new equilibrium points on the line: x — 1/2 — /i for Model A. If there is any 
new equilibrium points for this model, it must be on the x-axis. 



Let 



_ (1 -fi)(x + fi) /x(x + /x- 1) 
— ~ , r — — X, 



\X + fj,\ 



\x + jJ, — 1| 3 



(25) 
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and thus, 



P2(X) 



+ {x+^-i) 2 x > for a; > 1 fx, 
^$ ~ {x 4- 1? ~x, for - Li < x < 1 - Li, 

~ (x+zT) 2 ~~ (x+^-1) 2 ~~ X > for X < —Li. 



(26) 



We also set 



xM b 



(27) 



( x 2 +T 2)3/2- 



It is obvious that P2^x) contains the terms related to the gravitational forces from the 
central binary and the centrifugal force. Q2^x) is simply the term contributed by the belt. 
By defining these two functions, we have both the advantages that the physical meanings 
are clear and some properties about the new equilibrium can be proved easily. 

From Eq.(20), we have f(x, 0) = k(x) = —P2(x) + Q2(x). The equilibrium points (x e , 0) 
satisfy f(x e , 0) = k(x e ) = 0, that is, P2(x e ) = Q2(x e ). There are two properties for this part 
of our results: Property 3.4 is for the case of non-equal mass binaries, i.e. 0.5 > /i > 
and Property 3.5 is for the case of equal mass binaries, i.e. 11 = 0.5. In Property 3.4, we 
find that there are three equilibrium points, (xj,0) for i = 1,2,3, where x\ G (—00, —Li), 
x 2 G (0, 1 — Li) and x 3 G (1 — fi, 00). However, the results for the equilibrium points which 
x— coordinate is in the region (—Li, 0) are more complicated. As we can see in Property 3.4 
(C), there are two equilibrium points which x-coordinates are in the region (— /x,0), and in 
Property 3.4 (D), there is no equilibrium point with x— coordinate in the region (— fj.,0). 
We will see that, in fact, Property 3.5 (A) could be combined with Property 3.4 (B). 
However, Property 3.5 (B) studies the same condition as in Property 3.4 (C) but the results 
for that condition are completely different. Moreover, There is nothing like Property 3.4 
(A) and 3.4 (D) in Property 3.5. These differences are due to that P 2 ( _ ^) = — P 2 {x) and 
Q 2 {— x) = —Q2(x) when /i = 1/2 but the symmetry is broken when /i < 1/2. We therefore 
decided to separate Property 3.4 and 3.5. For convenience, we will denote /i + to represent 
that x tends to \x from the right hand side and /i~ to represent that x tends to ji from the 
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left hand side. 

Property 3.4 

(A) There is an x\ > 1 — fi and an x 2 G (0, 1 — ji) such that k(xi) = and k(x 2 ) = 0. 

(B) There is an x 3 < —fi such that k(x 3 ) = 0. 

(C) IfT < V2fi and g 2 (-T/ v / 2) > P 2 (-T/y/2), then there are two points in the region 
(— fi, 0) such that k(x) = 0. 

(D) 7/g 2 (-T/ v / 2) < P 2 (0), then there is no point in the region (—fi, 0) such that k(x) = 0. 
Proof of (A): From Eqs. (25)-(27), we know that if x > 1 — fi, then 



and 



Since P^x) < for a; > 1 — fi, P 2 is a monotonically decreasing function for 
x G (1 — fi, oo). Because lim x _>(i_ M )+ P 2 (x) = oo, linx^oo P 2 (rr) = — oo and Q 2 (l — fi) < 0, 
lim^ooQ^rr) = 0, we have lim x _ (1 _ Jll )+ = lim x _ > ( 1 _ /1 )+(-P 2 (x) + Q 2 (x)) < and 
liHia-^oo k(x) = lim 2 .^ 00 (— P 2 (x) + Q 2 (x)) > 0. Therefore, there is a point x 1 > 1 — fi such 
that fc(xi) = 0. 

In the region(0, 1 — fi), P 2 (0) > 0, lim ;i; _ > ( 1 _ /i )- P 2 (x) = — oo, (^2(0) = and 
Q 2 (l — fi) < 0, so we have k(0) = —P 2 (0) + Q 2 (0) < and lim a; ^( 1 _ jU )- = 00. Therefore, 
there is an x 2 G (0, 1 — fi) such that k(x 2 ) = 0. □ 

Proof of (B): If x < —fi, then 

(x + fi) 1 (x + fi — xy 

Since P 2 (x) < for x < —fi, P 2 is a monotonically decreasing function for x < —fi. 
Because lim x ^_ M - P 2 (x) = —00, lim 2 .^_ 00 P 2 (x) = 00 and Q 2 (—fi) > 0, 
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lim-r-^-oo Q2(x) = 0, we have linx^.^- k(x) > and \im x ^_ OQ k(x) < 0. Therefore, 
there is a point x 3 < —[/, such that k(x 3 ) = 0. □ 

Proof of (C): If x G {-Li, 0), then 

r> / \ 1 "A* A* 
r^lX) = — - — r— — X. 

V ; (x + /i) 2 (x + fi-1) 2 

If T < y/2fi, we have that —ii < —T/ y/2. We then discuss the possible points such that 
k(x) = for x = -T/y/2, x G {-T/y/2, 0) and x G (— //, —T/y/2) separately. 

If g 2 (-T/ v / 2) > P 2 (-T/ v / 2), then k(-T/y/2) = -P 2 (-T/y/2) + Q 2 {-T/y/2) > 0. 
Thus, fc(a;) ^ when x = -T/y/2. 

Because we set < \x < 0.5 in this paper, we have P 2 (0) > 0. Since Q 2 (0) = 0, we have 
fc(0) = -P 2 (0) + Q 2 (0) < 0. Further, if Q 2 (-T/ v / 2) > P 2 (-T/ v / 2), then k(-T/y/2) > 0. 
Thus, there is an x 4 G (— T/y/2, 0) such that fc(x4) = 0. 

On the other hand, since P£(x) < 0, Q' 2 (x) > for x G (-/I, -T/y/2) , 
k'(x) > for x G (— /i, -T/y/2). Thus, fc(a;) is a monotonically increasing function. 
Because lim. E ^_ M + P 2 (x) = 00 and Q 2 (-/■*) > 0, we have lim x _ > _ /1 + fc(a;) < 0. If 
g 2 (-T/ v / 2) > P 2 (-T/ v / 2), then k{-T/y/2) > 0. Therefore, there is a unique point 
x 5 G (-//, -T/v 7 ^) such that k(x 5 ) = 0. □ 

Proo/ 0/ (P)j: Since P 2 (x) < for x G (— /i,0), P 2 is a monotonically decreasing 
function and has a minimum value P 2 (0) in the region (— //, 0). It is easy to show that 
g 2 (-T/ v / 2) is the maximum value for Q 2 (x) in the region {—Li, 0). If Q2(—T/y/2) < P2(0), 
we have P 2 (x) > Q2(x) and thus k(x) < in the region (—fx, 0). Therefore, there is no point 
such that k(x) = in the region (—fx, 0). □ 

In Property 3.5, we study the fx — 1/2 case. Since P 2 (— x) = —P 2 (x) and 
Q 2 (— x) = —Q2(x), we only need to discuss (— 00, 0) region. 
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Property 3.5 

(A) There is an X\ < —1/2 such that k(xi) = 0. 

(B) IfT < l/y/2 and Q 2 {-T/y/2) > P 2 (-T/V2), then there is a point in the region 
(-1/2, -T/V2) such that k(x) = 0. 

Proof : (A) can be proved easily by the the same method as the proof for Property 
3.4 (B). In the following, we would like to prove (B). Since Q 2 (-T/V2) > P 2 (-T/y/2), 
k(—T/\/2) > 0. Since Q 2 (—l/2) is finite and lim x _>(_i/2)+ P 2 (x) = oo, lim x _ > (_ 1 / 2 )+ k(x) < 0. 
Thus, there is a point x* e (-1/2, -T/V2) such that k(x*) = 0. □ 

From Property 3.5, if T < 1/V2 and Q 2 (-T/^/2) > P 2 (-T/V2), we could have two 
points such that k(x) = in the region (— oo, 0). By the property of symmetry, we will 
have another two roots in the region (0, oo). Since k(0) = 0, there are at least five points 
such that k(x) = for the equal mass case. We know that there are only three equilibrium 
points on the x-axis for the classical restricted three-body problem, thus new equilibrium 
points exist for this case. From Property 3.3, 3.4 and Property 3.5, we have the following 
corollary to analytically determine the area which certainly have new equilibrium points for 
both < fx < 0.5 and /i = 0.5 cases. 

Corollary 3.6 

IfT < y/2fi and Q 2 {-T/y/2) > P 2 (-T/ v / 2), then we will have new equilibrium points. 
That is, there are at least five equilibrium points on the x-axis for System (10). 

One can see that the conditions studied in Property 3.4 and 3.5 do not include all 
possible cases. That is, there are some conditions we fail to provide analytical results due 
to that the results of those cases would more strongly depend on the detail structure of the 
function P 2 (x) and Q 2 (x). Thus, Corollary 3.6 shows that the condition, T < \^2fi and 
Q2(—T/V2) > P2(—T/y/2), is only a sufficient condition to have new equilibrium points. It 
is not a necessary condition. 
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To check the above statements, we numerically solve both Eq.(20) and Eq. (21) and 
find out the number of equilibrium points for different given parameters. The numerical 
scheme of root finding is the Van Wijngaarden-Dekker-Brent Method (Brent 1973). This is 
an excellent algorithm recommended by Press et al. (1992). We set a high level of accuracy 
that the maximum error is 10~ 8 for the locations of equilibrium points on both x-axis and 
the line: x — 1/2 — All the new equilibrium points are on the x-axis as we have proved 
in the analytic results. Figure 1 is the results on the \x — M b plane. That is, we numerically 
search the solutions of both Eq.(20) and Eq.(21) for different M b ), where \i is between 
[0,0.5] and M b is between [0, 0.6]. Those (//, M b ) with new equilibrium points are marked 
by points. 

In each panel, there is a solid curve which satisfies Q 2 {-T/y/2) = P 2 (-T/ v / 2) with 
given parameters. Figure 1(a) shows that this curve perfectly match the boundary of the 
existence area of new equilibrium points. Therefore, Corollary 3.6 precisely predicts our 
numerical results. However, as we mentioned earlier, our analytic results do not include all 
possible conditions. Thus, the solid curve and the boundary of the existence area of new 
equilibrium points do not match that well in Figure l(b)-(d). We cannot see the solid curve 
in Figure 1(d) since this curve is out of our studied parameter space. 

Moreover, Figure 1 shows that the area without any points are on the left-bottom of 
/i — M b plane. This is reasonable since (1) when M b is small, it becomes classical restricted 
three-body problem and (2) when /x is small, this system become similar to a planetary 
system with a belt. 

From Figure 1(a), i.e. T = 0.01 case, it is interesting that there is a lower limit of mass 
ratio: below that ratio, there would be no new equilibrium point no matter how large the 
belt mass is. Such limits of mass ratio become larger for larger T as in Figure (b)-(d). In 
general, when the mass of the belt is larger, it is easier to have new equilibrium points for 
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larger mass ratio. It is interesting that the new equilibrium points could exist for the larger 
area on the /x — M b plane when T is smaller. It is also interesting that the locations and 
number of equilibrium points depend on T but not directly on a and b. T is the summation 
of a (the flatness parameter) and b (the core parameter) and the Miyamoto-Nagai potential 
would be equivalent to the potential of a point mass when T = 0. 

Figure 2 is the results on the T — M b plane. That is, we numerically search the solutions 
of both Eq.(20) and Eq.(21) for different (T, M b ), where T is between [0,0.4] with grid size 

0. 01 and M b is between [0, 0.6] with grid size 0.02. Those (T, M b ) with new equilibrium 
points are marked by points. 

Figure 2 shows that when T is small enough, we always have new equilibrium points 
for any value of belt mass. We also plot out the solid curve: Q 2 (-T/ v / 2) = P 2 (-T/y/2) 
in each panel. The left-up side of the solid curve satisfies the condition of Corollary 3.6, 

1. e. Q 2 {— T/y/2) > P 2 { — T/y/2), so we have new equilibrium points in this area. We find 
that these solid curves can match the boundaries of the existence area of new equilibrium 
points very well for Figure 2(a)-(c) but not that well for Figure 2(d). Nevertheless, they all 
indicate that our analytic results are completely consistent with our numerical results. 

In addition to the full circle points, for the comparison purpose here, there are also open 



triangle points in Figure 2(d). The open triangle points are the results for n = y 1 — 2/ 6 (0.5) 
as in Jiang & Yeh (2006). The open triangle points almost cover the same area as the one 
covered by the full circle points, thus our simplification that n = 1 is a good approximation. 




3.2. 



Model B 



On the other hand, we also study Model B: Power-Law Profile. By Property 3.1, the 
equilibrium points (x e , y e ) must be on the lines: x — 1/2 — fi or y = 0. Property 3.1 and 3.2 
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indicates that we only need to solve x e from 



x 



k(x) 
and y e from 
h{y) = 1 



1 



Zx f r ° 

— / p(r')r' 



E F 
+ 



ny ryf ny j ryf 



dr' 



y=0 



(29) 



-3/2 2 

— / p(r')r' 



r 2 j i 



E 



+ 



ny ryf ny j ryf 



dr' 



= 



X=(/li-/12)/2 



(30) 



to get all the equilibrium points (x e , 0) and (1/2 — p, y e ). Because there are Elliptic Integrals 
in k(x) and h(y), it is very difficult to get any analytic results as Property 3.3, 3.4 and 
3.5 of Model A. Thus, we try to find the roots of k(x) and h(y) numerically and we only 
present the numerical results here. The numerical method is exactly the same as the one 
used in Model A. 

Contrary to Model A, we find that the new equilibrium points can exist both on the 
vertical line: x — 1/2 — p: and on the rr-axis. Figure 3 and Figure 4 are the results, where 
the cross points are those area when there are more than three roots for k(x) = 0, the circles 
are the area when there are more than two roots for h(y) = 0. Please note that for the 
classical restricted three-body problem, there are three collinear points and two triangular 
ones. 

Figure 3 is the results on the fx — M b plane. That is, we numerically search the solutions 
of both Eq.(20) and Eq.(21) for different (fi,M b ), where fi is between [0,0.5] and M b is 
between [0, 0.6]. 

The area of circles in Figure 3(a) looks like a square, thus there is a lower limit of mass 
ratio that no matter how large the belt mass is, there would be no new equilibrium point. 
There is also a lower limit of belt mass that no matter what mass ratio is, there would be 
no new equilibrium point. Figure 3(b) and (c) are very similar but the area becomes more 
like a quarter of a circle. The existence area of new equilibrium points in Figure 3(d) is 
very different from others. 
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The cross points in Figure 3 only appear for a small range of mass ratio. That is, with 
the additional force from the belt, the mass ratio has to be a particular value to have new 
equilibrium points. This means that it is more difficult to have new equilibrium points 
on the x-axis than that on the y-axis. To have new equilibrium points on the y-axis, the 
centrifugal force has to balance with the y-component of the total force from both central 
binary objects plus the force from the belt. That is, there are three components to balance 
each other at the equilibrium points and the mass ratio would only modify the locations of 
equilibrium points. However, on the x-axis, the forces from two central binary objects are 
in different directions at new equilibrium points, thus there are four components to balance 
each other for this case. Therefore, it will be more difficult for them to get balanced. To 
have new equilibrium points, the mass ratio has to be within a particular range. 

Figure 4 is the results on the — M b plane. That is, we numerically search the 
solutions of both Eq.(20) and Eq.(21) for different (r i: M b ), where r { is between [0,1] with 
grid size 0.01 and M b is between [0, 0.6] with grid size 0.02. 

The cross points in Figure 4 show that there is a gap of r\ where it is more difficult 
to have new equilibrium points on the x-axis. It seems that this gap become larger when 
/i is larger. The reason why it would be more difficult to have new equilibrium points for 
particular region of r« might be related to the relative positions of central binary objects 
with respect to the belt. The number of separate components which could balance each 
other might influence the existence of new equilibrium points. 

The circles in Figure 4 are very interesting. In Figure 4(a) where the mass ratio is 
fx = 1/11, there is a triangle area, i.e. when 0.5 < r\ < 0.8, the range of belt mass to have 
new equilibrium points is larger if is larger. However, when r« < 0.5 or r« > 0.8, there 
is no new equilibrium point. The triangle becomes slightly larger in Figure 4(b) and it 
becomes the combination of two trapezoids facing each other in Figures 4(c)-(d). 
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4. Stability Analysis 



When the locations of equilibrium points are determined, it would be interesting to 
understand the stability properties around these points. We now consider the following 
system: 



dx _ 
dt ~ u 

dt u 



% = 2v + f(x,y) 



§ = -2u + g(x,y), 



(31) 



where f(x,y), g(x,y) are defined in Eqs.(15)-(16). Following the usual linearization, the 
eigenvalues A corresponding to the equilibrium points would satisfy: 



A 4 + (4 - f x - g y )X 2 + 2(f y - g x )X + f x g y - g x f y = 0, 
where f x = df(x,y)/dx and f y ,g x ,g y are also defined similarly. 



(32) 



Consider a new equilibrium point (x e ,y e ) of Model A as an example, so (x e ,y e ) satisfies 
k(x e ) = and y e = 0. Thus, we have: 

3M b x 2 e 



g y (x e ,0) = 



Ob g J ^Qljq I jjj I 

nii-n)( i 



|x e + /i-l| 3 (x 2 e + T 2 ) 5 / 2 ' 



1 



\X e — /i| 3 \X e + jJL — 1| 



(33) 



Since f y (x e ,0) = and g x (x e: 0) = 0, Eq.(32) becomes 

\ 4 +(4-f x -g y )\ 2 + f x g y = 0. 
For convenience, we define Q = f x g y and U — f x + g y — 4. Therefore, we have 



(34) 



x9 n + vn 2 - 4Q , x9 n - vn 2 - 4^ 

X + = and A - 



2 - 2 

These two relations will determine the eigenvalues and the properties of the equilibrium 
point. In fact, the signs of II 2 - 4fi, II + y/Il 2 - 4Q, and II - VII 2 - 4fi will completely 



(35) 
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determine the results. However, when one knows the signs of II and fi, the details of possible 
eigenvalues can be directly determined by them. We therefore list all the combinations in 
Table 1, in terms of the signs of II 2 — 4fl, IT and Q. 



Table 1. The Possible Combinations 



Condition A 


Condition B 


Condition C 


Row 




A 2 


n 2 - 40 < 






1 


c 


c 


n 2 - 40 = o 


n < 




2 


— 


— 


n 2 - 40 = o 


n = o 




3 








n 2 - 40 = o 


n > o 




4 


+ 


+ 


n 2 - 40 > o 


n < o 


< 


5 


+ 




n 2 - 40 > o 


n < o 


= 


6 







n 2 - 40 > o 


n < o 


> 


7 






n 2 - 40 > o 


n = o 


< 


8 


+ 




n 2 - 40 > o 


n = o 


= 


9 (impossible) 






n 2 - 40 > o 


n = o 


> 


10 (impossible) 






n 2 - 40 > o 


n > o 


< 


11 


+ 




n 2 - 40 > o 


n > o 


= 


12 


+ 





n 2 - 40 > o 


n > o 


> 


13 


+ 


+ 



In Row 1, due to II 2 — 4f2 < 0, A^ are then imaginary complex conjugate numbers. 
Thus, A± can be expressed as A + = ±a±bi and A_ = ±c±di, where a,b,c and d are positive 
real numbers. So in this case, the equilibrium point is an unstable point. 

In Row 2, 3, and 4, 

Thus, A± could be pure imaginary (Row 2), zeroes (Row 3) and ±^Af (Row 4). The 
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equilibrium point is neutrally stable for Row 2 and 3 but unstable for Row 4. 

In Row 5, because tt < 0, we have IT 2 - 4fi > n 2 . Thus, Vn 2 - 4fi > | XT | = -II. In 
this case, II + Vn 2 - 4fi > and II - ^Yi 2 - AVt < 0. We then know that A + = ±a, where a 
is a positive real number; A_ = ±di, where d is a positive real number. So, the equilibrium 
point is an unstable point. 

In Row 6, because f2 = 0, we have 

n± mi 



A 2 = 



2 

Because IT < 0, A+ = and X 2 _ — IT < 0. Thus, A + = and A_ is pure imaginary. So in 
this case, the equilibrium point is a neutrally stable point. 

In Row 7, because Q > 0, we have II 2 -4Q < II 2 . Therefore, Vn 2 - 4f2 < | XT | = -II. In 
this case, II + Vn 2 — AVt < and II — \/n 2 — 4f2 < 0. Both A + and A_ are pure imaginary. 
So, the equilibrium point is a neutrally stable point. 

In Row 8, because II = 0, we have 



K = ^ and ,> = z^. 

+ 2 "2 
Since Q < 0, thus A^ > and A^ < 0. We then have that A + = ±a, where a is a positive 
real number; A„ = ±di, where d is a positive real number. So, the equilibrium point is an 
unstable point. 

In Row 9 and 10, Condition C contradicts with Condition A after we use Condition B: 
II = 0. Thus, these two are impossible. 



In Row 11, because tt < 0, we have IT 2 - 4fi > II 2 . Thus, VU 2 - AVI > |IT| = IT. In this 
case, II + Vn 2 - 4fi > and II - Vn 2 - 4fi < 0. We then know that A + = ±a, where a is 
a positive real number; A_ = ±di, where d is a positive real number. So, the equilibrium 
point is an unstable point. 
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In Row 12, because f2 = 0, we have 

= n±jn|. 

± 2 

Because IT > 0, A5j_ = II and A 2 = 0. Thus, A+ = iv 7 !! and A_ = 0. So in this case, the 
equilibrium point is unstable. 

In Row 13, because Q > 0, we have IT 2 - 4fi < II 2 . Therefore, Vn 2 - 4fi < |IT| = n. 
In this case, II + Vn 2 — 4f2 > and II — Vn 2 — 4f2 > 0. Both A+ and A_ are real and one 
of their roots is positive. So, the equilibrium point is unstable. 

In general, from the above analysis, there are two categories: (i) when all the four 
eigenvalues' real parts are zero, the equilibrium point is neutrally stable, (ii) as long as 
one of the eigenvalues has a positive real part, the equilibrium point is unstable. For this 
system, it is impossible that all the four eigenvalues' real parts are negative. In Table 1, 
Row 2, 3, 6, 7 belong to the first category because their eigenvalues are either zero or pure 
imaginary. Row 1, 4, 5, 8, 11, 12, 13 belong to the second category because at least one of 
the eigenvalues is positive. 

To make it convenient to use Table 1 in the future, the properties of A 2 ^ and A 2 are 
indicated in the final two columns of the table. In that two columns, "c" means imaginary 
complex conjugate numbers, "+" means positive numbers, "— " means negative numbers, 
and "0" means zeroes. 

As an example, when we set M b = T = 0.01 and /i = 4/9 in Model A, there are 
five equilibrium points on rr-axis: (-1.180392,0), (-0.060183,0), (-0.000137,0), (0.118920,0), 
(1.218591,0). We find that only (-0.000137,0) belongs to the first category and is neutrally 
stable. All the others belong to the second category and is unstable. 
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5. Concluding Remarks 

We have studied the existence of new equilibrium points for our system. In addition 
to the usual equilibrium points, we find that there are some new equilibrium points for 
particular parameter space. We provide two models in this paper. One is an analytic model 
from Miyamoto & Nagai (1975) (Model A) and another is a model of power-law density 
profile (Model B). Since there is no analytic formula for the power-law density profile, we 
have to study this model numerically. For Model A, we have some interesting analytic 
results and they are consistent with the numerical results. The new equilibrium points only 
exist on the rr-axis for Model A. However, for Model B, we can have new equilibrium points 
both on the vertical line: x — 1/2 — /i and the x-axis. 

Our study implies that the structure of constant potential contour might be more 
complicated due to the existence of new equilibrium points. The stability could be different 
according to the linear analysis in the previous section. Thus, the dynamics around a binary 
system might become very different when the influence of a belt is included. 
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Fig. 1. — The existence area of new equilibrium points on the \i — M& plane for Model A. In 
(a) T = 0.01; (b) T = 0.1; (c)T = 0.2 and (d) T = 0.3 (please see the main text for more 
details). 
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Fig. 2. — The existence area of new equilibrium points on the T — M& plane for Model A. In 
(a) ji = 1/11; (b) \x = 1/6; (c)/x = 1/3 and (d) // = 1/2 (please see the main text for more 
details). 
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Fig. 3. — The existence area of new equilibrium points on the [i — Mb plane for Model B. 
In (a) Ti = 0.2; (b) n = 0.4; (c)rj = 0.6 and (d) r\ = 0.8 (please see the main text for more 
details). 
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Fig. 4. — The existence area of new equilibrium points on the — M& plane for Model B. In 
(a) /j, = 1/11; (b) \x = 1/6; (c)/i = 1/3 and (d) /i = 1/2 (please see the main text for more 
details). 



